Confined sandpile in two dimensions: Percolation and singular diffusion.

We investigate the properties of a two-state sandpile model subjected to a confining potential in two dimensions. From the microdynamical description, we derive a diffusion equation, and find a stationary solution for the case of a parabolic confining potential. By studying the systems at different confining conditions, we observe two scale-invariant regimes. At a given confining potential strength, the cluster size distribution takes the form of a power law. This regime corresponds to the situation in which the density at the center of the system approaches the critical percolation threshold. The analysis of the fractal dimension of the largest cluster frontier provides evidence that this regime is reminiscent of gradient percolation. By increasing further the confining potential, most of the particles coalesce in a giant cluster, and we observe a regime where the jump size distribution takes the form of a power law. The onset of this second regime is signaled by a maximum in the fluctuation of energy.

Fragmentation processes in two-phase materials.

We investigate the fragmentation process of solid materials with crystalline and amorphous phases using the the discrete element method. Damage initiates inside spherical samples above the contact zone in a region where the circumferential stress field is tensile. Cracks initiated in this region grow to form meridional planes. If the collision energy exceeds a critical value which depends on the material's internal structure, cracks reach the sample surface resulting in fragmentation. We show that this primary fragmentation mechanism is very robust with respect to the internal structure of the material. For all configurations, a sharp transition from the damage to the fragmentation regime is observed, with smaller critical collision energies for crystalline samples. The mass distribution of the fragments follows a power law for small fragments with an exponent that is characteristic for the branching merging process of unstable cracks. Moreover this exponent depends only on the dimensionally of the system and not on the microstructure.

Scaling laws for impact fragmentation of spherical solids.

We investigate the impact fragmentation of spherical solid bodies made of heterogeneous brittle materials by means of a discrete element model. Computer simulations are carried out for four different system sizes varying the impact velocity in a broad range. We perform a finite size scaling analysis to determine the critical exponents of the damage-fragmentation phase transition and deduce scaling relations in terms of radius R and impact velocity v(0). The scaling analysis demonstrates that the exponent of the power law distributed fragment mass does not depend on the impact velocity; the apparent change of the exponent predicted by recent simulations can be attributed to the shifting cutoff and to the existence of unbreakable discrete units. Our calculations reveal that the characteristic time scale of the breakup process has a power law dependence on the impact speed and on the distance from the critical speed in the damaged and fragmented states, respectively. The total amount of damage is found to have a similar behavior, which is substantially different from the logarithmic dependence on the impact velocity observed in two dimensions.

Fragmentation processes in impact of spheres.

We study the brittle fragmentation of spheres by using a three-dimensional discrete element model. Large scale computer simulations are performed with a model that consists of agglomerates of many particles, interconnected by beam-truss elements. We focus on the detailed development of the fragmentation process and study several fragmentation mechanisms. The evolution of meridional cracks is studied in detail. These cracks are found to initiate in the inside of the specimen with quasiperiodic angular distribution. The fragments that are formed when these cracks penetrate the specimen surface give a broad peak in the fragment mass distribution for large fragments that can be fitted by a two-parameter Weibull distribution. This mechanism can only be observed in three-dimensional models or experiments. The results prove to be independent of the degree of disorder in the model. Our results significantly improve the understanding of the fragmentation process for impact fracture since besides reproducing the experimental observations of fragment shapes, impact energy dependence, and mass distribution, we also have full access to the failure conditions and evolution.

Universality behind Basquin's Law of Fatigue.

Basquin's law of fatigue states that the lifetime of the system has a power-law dependence on the external load amplitude, tf approximately sigma 0- alpha, where the exponent alpha has a strong material dependence. We show that in spite of the broad scatter of the exponent alpha, the fatigue fracture of heterogeneous materials exhibits universal features. We propose a generic scaling form for the macroscopic deformation and show that at the fatigue limit the system undergoes a continuous phase transition. On the microlevel, the fatigue fracture proceeds in bursts characterized by universal power-law distributions. We demonstrate that the system dependent details are contained in Basquin's exponent for time to failure, and once this is taken into account, remaining features of failure are universal.

Computer simulation of fatigue under diametrical compression.

We study the fatigue fracture of disordered materials by means of computer simulations of a discrete element model. We extend a two-dimensional fracture model to capture the microscopic mechanisms relevant for fatigue and we simulate the diametric compression of a disc shape specimen under a constant external force. The model allows us to follow the development of the fracture process on the macrolevel and microlevel varying the relative influence of the mechanisms of damage accumulation over the load history and healing of microcracks. As a specific example we consider recent experimental results on the fatigue fracture of asphalt. Our numerical simulations show that for intermediate applied loads the lifetime of the specimen presents a power law behavior. Under the effect of healing, more prominent for small loads compared to the tensile strength of the material, the lifetime of the sample increases and a fatigue limit emerges below which no macroscopic failure occurs. The numerical results are in a good qualitative agreement with the experimental findings.

Counterion correlations and attraction between like-charged macromolecules.

A simple model is presented for the appearance of attraction between two like-charged polyions inside a polyelectrolyte solution. The polyions are modeled as rigid cylinders in a continuum dielectric solvent. The strong electrostatic interaction between the polyions and the counterions results in counterion condensation. If the two polyions are sufficiently close to each other their layers of condensed counterions can become correlated resulting in attraction between the macromolecules. To explore the counterion induced attraction we calculate the correlation functions for the condensed counterions. It is found that the correlations are of very short range. For the parameters specific to the double stranded DNA, the correlations and the attraction appear only when the surface-to-surface separation is less than 7 A.